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We present a detailed analysis of the distance determination for one Intermediate Velocity 
Cloud (IVC G139.6 + 47.6) from the ongoing DDO IVC Distance Project. Stars along the line 
of sight to G139.6 + 47.6 are examined for the presence of sodium absorption attributable to the 
cloud, and the distance bracket is established by astrometric and spectroscopic parallax mea- 
surements of demonstrated foreground and background stars. We detail our strategy regarding 
target selection, observational setup, and analysis of the data, including a discussion of wave- 
length calibration and sky subtraction uncertainties. We find a distance estimate of 129 ± 10 
pc for the lower limit and 257I33 1 f° r ^ ne upper limit. Given the high number of stars showing 
absorption due to this IVC, we also discuss the small-scale covering factor of the cloud and the 
likely significance of non-detections for subsequent observations of this and other similar IVC's. 
Distance measurements of the remaining targets in the DDO IVC project will be detailed in a 
companion paper. 

Subject headings: ISM: clouds — ISM: individual (G139.6 + 47.6) — ISM: structure — stars: distances 
— techniques: spectroscopic 

1. INTRODUCTION behind the cloud, one hopes that metallic lines at 



ABSTRACT 



As with many objects in the Universe, Interme- 
diate Velocity Clouds (IVC's) are poorly under- 
stood due to a lack of distance estimates. Many 
of the interesting physical properties of the IVC's 
(mass, density, volume, etc.) depend on the dis- 
tance, and without these properties it is difficult 
to constrain models concerning their origins, evo- 
lution, or relative importance to the interstellar 
medium (Wakker & van Woerden 1997). 



the cloud's systemic velocity can be detected in the 
star's spectrum. By determining whether stars are 
behind the cloud, we set upper and lower limits on 
its distance. Astrometric or spectroscopic paral- 
lax can then be used to determine the distances 
to the observed stars to bracket the cloud's dis- 
tance. Such a technique, unfortunately, requires a 
large investment of telescope time (Wakker & van 
Woerden 1997). 



There are several ways to determine the dis- 
tances to IVC's, which have been used to vary- 
ing degrees of success (see Gladders et al. 1998, 
and references therein). One of these methods 
stands out as particularly successful in determin- 
ing both upper and lower limits on the distances 
to the clouds: the absorption line method. In 
this method, stars along the line of sight to the 
cloud are used as distance markers. If a star is 



We began the DDO IVC Distance Project in 
the summer of 1997 in the hopes that a reliable 
distance to the Draco cloud could be determined 
(see Gladders et al. 1998). The success of this 
project soon led to distance determinations for 
other IVCs (see Clarke et al. 1999). In this paper, 
we use the cloud G139.6 + 47.6 as a case study in 
order to present our methodology as a precursor to 
the publication of the entire data set, forthcoming 
in a companion paper. This cloud's Local Stan- 
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dard of Rest velocity is t>LSR = — 12 km s _1 and 
is one of the lowest in our sample, which range 
from 3.5 km s _1 to —47 km s _1 . Traditionally, 
IVC's are distinguished by wlsrI < 90 km s _1 
and G139.6 + 47.6, having a relatively low veloc- 
ity, might be more appropriately labeled a Low 
Velocity Cloud (LVC). However, as we are not dis- 
cussing the origins or kinematics of these clouds, 
we will not make this distinction here. 

2. OBSERVATIONS 
2.1. Target Selection 

We have taken our target list of clouds from the 
sample of Heiles, Reach, & Koo (1988). Among 
the 26 clouds, we selected 16 to observe that are 
at high galactic latitudes and are observable from 
the David Dunlap Observatory. Several of these 
clouds have been grouped together in pairs, as 
they do not seem to be separated either in an- 
gular or velocity space. We are thus left with a 
sample of 11 cloud complexes. 

With the clouds selected, we then require good 
stellar targets along their lines of sight. We con- 
sider a "good" target to be a star with B — V < 
0.35, as we want early-type stars with as few stel- 
lar features as possible so that absorption features 
can be easily identified as originating in the cloud 
rather than the stellar photosphere. We use the 
Tycho Catalogue (Hoeg et al. 1997) and the USNO 
Catalog (Urban et al. 1998) to select our targets. 
Cloud sizes are estimated through both the H I 
atlas of Hartmann & Burton (1997) and through 
IRAS 100 fi maps (see Figure 1). We then search 
the Tycho Catalogue, using the online tools, for 
stars with positions within the cloud and with 
the appropriate colors. When the positions of the 
stars have been determined, we build an H I col- 
umn density profile in velocity space at the po- 
sition of each star (see Figure 2). Assuming hy- 
drogen or IRAS emission is a good indication of 
sodium absorption (the validity of this will be dis- 
cussed in §4.2), we can then determine which stars 
represent good probes of the cloud. Also, using 
the B — V color index, the stars' distances can be 
estimated, which further helps in target selection. 
Figure 1 shows the IRAS 100 fi emission for the 
IVC cloud G139.6 + 47.6 as well as the positions of 
the 16 target stars that satisfied the above criteria. 



2.2. Observational Setup 

Since the beginning of this project, we have 
been assigned a total of 137 nights, of which 63 
have produced useful spectra for at least a por- 
tion of the night. To date we have observed 86 
stars along lines of sight to eleven different IVC's. 

Nominally, each target star is observed using 
two observational setups: a sodium setup and a 
classification setup. Both use the David Dun- 
lap Observatory (DDO) 1.88 m telescope and 
Cassegrain spectrograph. The classification setup 
is used to classify each star so that its distance 
can be determined. The sodium setup is used to 
observe the Na I doublet to determine whether the 
star is behind the cloud or not. The sodium setup 
is at higher resolution than the classification setup, 
and so it represents the bulk of the observations. 
In addition to our program stars, spectroscopic 
standards are observed in both the sodium and 
classification setups. We observe standards in the 
sodium setup in order to generate template spec- 
tra for use with the cross-correlation techniques, 
which are detailed in §4.1. 

2.2.1. Sodium Setup 

The sodium setup yields a wavelength coverage 
of approximately 5800-6000 A with a resolution 
of 0.43 A. With this setup, we cannot resolve any 
sodium absorption caused by the cloud (our veloc- 
ity resolution is 22 km s _1 at the Na I doublet); 
however, this is not necessary for our purposes (see 
§4.1). We observe targets down to a magnitude of 
12.5, for which it takes 4 hours to achieve a sig- 
nal to noise ratio of 50. The wavelength region of 
interest is significantly affected by water vapor ab- 
sorption lines, so the observation of telluric stan- 
dards is important. Typically, we observe three 
to four telluric standards throughout the night. 
These standards are later used to subtract the tel- 
luric lines from the sodium spectra. 

2.2.2. Classification Setup 

For the classification setup, we cover the range 
3800-4400 A at a resolution of 2.5 A. We obtain 
classification spectra of any stars which show defi- 
nite absorption or lack of absorption by the cloud 
and that do not have Hipparcos distance measure- 
ments. Classification spectra of standards taken 
from Garcia (1989) are also obtained, covering 
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spectral types BO to G5 and luminosity classes I 
to V. These standards are used to determine the 
distances to several program stars using spectro- 
scopic parallax. 

3. DATA ANALYSIS 

The bulk of the data reduction is now done 
through a processing pipeline developed using 
standard IRAF tasks. This pipeline allows us 
to quickly reduce a night's data and also ensures 
a consistent procedure for all data. The rapid 
turnaround allows greater precision in subsequent 
target selection. Further, with over 20 observers 
working on the same project, a centralized repos- 
itory of data is necessary. Such a database, im- 
plemented in HTML, was devised and became an 
invaluable resource. 

The data are de-biased and flat-fielded using 
standard methods. Each two-dimensional spec- 
trum is then rectified to a common wavelength 
scale. With the DDO's close proximity to one of 
North America's largest urban centers, one must 
be very careful in the subtraction of the sodium 
sky lines. Even though most of the light from 
high-pressure sodium lights is self-absorbed, the 
two sodium lines from low pressure lamps are quite 
bright (see Figure 3). Great care is taken to ensure 
proper rectification of the 2-D spectra and hence 
proper sky line subtraction. 

Lastly, we remove the telluric lines using the 
spectrum of the telluric standard that is observed 
the closest in time to the program star. This 
is done by normalizing the reduced, continuum- 
subtracted 1-D telluric spectrum and then cor- 
recting for the difference in airmass. Figure 4 
shows a typical spectrum of the telluric stan- 
dard HD177724. Two telluric standard stars 
were chosen due to their lack of contaminating 
inter-stellar features in the wavelength region of 
interest: HD177724 (AOVn, m v = 2.988) and 
HD120315 (B3V, m v = 1.852). In order to sam- 
ple the slit effectively without saturating the CCD, 
neutral density filters were used. 

In order to test for the wavelength calibra- 
tion stability, the sky lines are extracted from 
the dispersion-corrected 2-D spectra of every pro- 
gram star. Using the same cross-correlation tech- 
niques outlined in §4.1, the shift in the sky Na I 
emission lines is computed. The average shift is 



0.35 km s _1 with a standard deviation of approx- 
imately 1 km s _1 , which is consistent with no sys- 
tematic shift and demonstrates the stability of the 
wavelength calibration. 

4. DISTANCE DETERMINATION 

There arc three steps to the determination of a 
cloud's distance: 1) testing for non-stellar contri- 
bution in Na I absorption lines, 2) assigning each 
star to the background or foreground, and 3) com- 
puting the distance to each star. Each of these 
steps is described in detail below. 

4.1. Detection of Absorption by the Cloud 

The most important step in the determina- 
tion of a cloud distance is the detection (or non- 
detection) of Na I from the cloud in the spectrum 
of stars along the line of sight. Ideally, one would 
like to use early-type stars to minimize the number 
of stellar lines in the vicinity of the Na I doublet, 
and one would like to resolve the absorption due 
to the cloud. However, to maximize the number of 
target stars for this project, we have chosen later- 
type stars and lower resolution than other authors 
(Benjamin et al. 1996). Because moderate disper- 
sion allows for a larger sampling of the underly- 
ing stellar spectrum, we can use cross-correlation 
techniques to decouple absorption due to the stel- 
lar envelope and any possible absorption due to 
the IVC. The technique is as follows. 

For each star in our sample, we estimate the 
spectral class using the sodium part of the spec- 
trum, or the classification spectrum if this has 
been obtained already. Based on this classifica- 
tion, we build a template sodium spectrum by 
measuring the central wavelength and equivalent 
width of all significant lines in the sodium spec- 
trum of a standard star with similar spectral 
class. The template spectrum is composed of a 
flat continuum with absorption lines at the same 
wavelengths and with the same equivalent widths 
(though with FWHM less than the resolution of 
the real spectra) as would be measured at rest 
(v r = 0). Using the IRAF task rvcorrect, each 
program star's spectrum is corrected to remove 
the velocity of the Earth and of the Sun relative to 
the Local Standard of Rest (LSR). Measurement 
of the Dopplcr shift in the absorption features rel- 
ative to their rest wavelengths will therefore yield 
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the LSR velocity of the star. 

Once the template spectrum has been pro- 
duced, we cross-correlate it with the spectrum of 
the program star in two separate regions. The first 
region, which we shall term the stellar region, is 
the whole sodium spectrum but excluding the re- 
gion around the Na I doublet. The second region 
we shall term the sodium region, which contains 
only the sodium doublet. The cross-correlation 
of the stellar region with the template spectrum 
gives us an estimate of the radial velocity and the 
average width of the stellar features. The cross- 
correlation of the sodium region with the template 
spectrum gives us an estimate of the velocity and 
average width of both the stellar features and the 
interstellar features if both are present in the re- 
gion. By comparing the cross-correlations, one can 
determine whether or not there are any Na I fea- 
tures due to interstellar absorption and, therefore, 
whether the star is behind the cloud. If the stel- 
lar features are too narrow and too close to the 
velocity of the IVC, the star is rejected. Figure 
5 shows the cross-correlations for the 4 program 
stars we have identified as background stars. The 
left panels show the Na I region of the spectra and 
the right panels show the cross-correlations with 
the template spectra. 

We consider a star to have a sodium detection 
if the cross-correlation of the Na I D lines peaks 
near the velocity of the cloud and either there are 
no significant stellar lines or any stellar lines are 
at a significantly different velocity than that of the 
cloud. There are cases where the velocity of the 
stellar lines is too close to the velocity of the cloud, 
yet the cross-correlation peak of the stellar lines 
is significantly broader than that of the sodium 
lines, which is suggestive of absorption due to the 
cloud. In these cases, we label the stars as pos- 
sible detections. Likewise, there are cases where 
the sodium and stellar lines have identical cross- 
correlation profiles near the velocity of the cloud 
and we label these as possible non-detections (see 
Table 1). 

4.2. Detection vs. Non-detection 

While the detection of absorption features due 
to the cloud in the spectrum of a star is strong 
evidence that the star is behind the cloud, the 
reverse is not necessarily true. A non-detection 
for a star can either mean the star is in the fore- 



ground, or it can simply mean the star is behind 
the cloud but the column density of Na I is not suf- 
ficient to produce detectable features in its spec- 
trum. Two uncertainties arise at this point. The 
first is the covering factor of the cloud on small an- 
gular scales. The H I column density is sampled 
from a relatively large beam, whereas our stars 
sample sodium absorption on an extremely small 
angular scale. It is therefore possible that small- 
scale variations in the covering factor are unre- 
solved by the large H I beam. The second uncer- 
tainty is the actual abundance of Na I relative to 
H I. In the absence of any Na I absorption, the ra- 
tio of Na I/H I is uncertain and could be very 
low. Therefore, clouds with no detections may 
very well have low Na I/H I. This problem could be 
resolved by extrapolation from other cloud metal- 
licities, once they are determined to some degree 
of accuracy. An example of low Na I/H I could be 
G86.5 + 59.6, which has eluded all our attempts 
at detecting Na I absorption. Our current lower 
limit to the distance to the cloud is 430 pc (Clarke 
ct al. 1999). However, its rather large radial ve- 
locity of —39 km s -1 may imply a large distance 
based on the terminal drag model of Benjamin 
et al. (1996), which predicts a linear correlation 
between the velocity perpendicular to the galactic 
plane and height above the plane. 

Recent measurements of IVC sight-lines have 
shown that Na I/H I varies on the order of 70 per 
cent over 30 arcsec and by as much as a factor 
of 2 over 1 arcmin (Smoker et al. 2002). Given 
the large angular distances between our target 
stars, we cannot expect the Na I/H I ratios ob- 
tained from our detections to be constant, nor can 
we use them to accurately predict the expected 
equivalent widths of Na I in other targets. This 
makes the identification of foreground stars more 
uncertain. Nevertheless, when evaluating the sig- 
nificance of a non-detection, the lowest observed 
EW(Na I) /Nhi for the cloud is used to predict the 
sodium absorption one would expect if the star 
were behind the cloud. If this absorption is de- 
tectable by our instruments yet is not observed, 
we conclude the star is in the foreground. Fig- 
ure 6 shows the expected sodium absorption in 
the spectrum of TYC 4151 1452 if the star were 
behind the cloud. The 4 spectra were computed 
using the EW(Na I) /Nhi obtained from our 4 de- 
tections. 
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To address the issue of covering factor, we ex- 
amine the distance distribution of our target stars. 
Of the 16 program stars, 4 have definite detections 
of Na I absorption. By obtaining distances to all 
these stars, one can now start to answer the ques- 
tion of covering factor, at least for this particular 
cloud complex. We wish to determine the prob- 
ability that a background star will show sodium 
absorption due to the cloud given the H I column 
density along its line of sight. This can be deter- 
mined in the following way. We first determine 
the best upper limit on the distance to the cloud 
(i.e., the closest star with a detection). We then 
look for any stars which could be farther away and 
count the number of detections versus the number 
of non-detections. If we find a non-detection which 
is farther away than a detection, this implies a cov- 
ering factor less than unity 

If we leave out the stars for which the detections 
and non-detections are of low confidence, then 
there are 4 stars with detections behind the nom- 
inal position of the cloud and one non-detection 
in front. If we include the possible detections and 
non-detections, then there are 7 detections behind 
the furthest non-detection and two non-detections 
in front of the closest detection. In both cases, 
this is consistent with a covering factor of unity. 

4.3. Stellar Distance Determination 

For the determination of the stellar distances, 
we use either trigonometric parallax from Hippar- 
cos or spectroscopic parallax using spectroscopic 
classification obtained either in the literature or 
from our own observations. Clearly, the faintest 
and farthest stars require spectroscopic classifica- 
tions to obtain accurate distances. However, since 
the goal is to determine the best distance to any 
given H I cloud, we are more concerned with accu- 
rately determining the distance to those stars that 
are closest to the cloud. 

To determine distances based on the spectral 
type of the star, we use the the spectral class and 
absolute magnitude calibration of Corbally & Gar- 
rison (1984) to assign an absolute magnitude to 
each classified star. The Tycho apparent magni- 
tude is then used to compute the distance mod- 
ulus. No correction for interstellar absorption is 
made as it is negligible at this galactic latitude and 
longitude, being near the Lockman Hole. We also 
ignore any extinction that may be caused by the 



cloud itself, as this is likely small (Stark, Kalberla, 
& Guesten 1997) and can only reduce the distance 
of the background stars, thereby improving our 
distance estimates. In assigning uncertainties to 
our distance estimates, two factors contribute the 
most: the uncertainty in the classification of the 
spectra and the intrinsic scatter in the absolute 
magnitude versus spectral type relation. For the 
former, we assume an uncertainty of plus or mi- 
nus one spectral class. For the latter, we take 
±0.7 mag (Jaschek & Gomez 1996). For the F- 
type stars, there is also a larger uncertainty in es- 
tablishing the luminosity class, which translates to 
a very large distance uncertainty. In these cases, 
we give distance ranges rather than uncertainties. 

5. RESULTS AND CONCLUSIONS 

Of the original 20 program stars, 4 were dis- 
carded due to strong stellar features at the veloc- 
ity of the cloud. Seven stars were inconclusive and 
require higher SNR and/or higher spectral resolu- 
tion than DDO can provide in order to make use of 
them. Table 1 lists the remaining 9 stars by Tycho 
name, the apparent magnitude, spectral type, dis- 
tance, equivalent width, H I column density due 
to the cloud, and note regarding detection. 

The DDO IVC Distance Project has been a re- 
sounding success and has far exceeded our initial 
expectations. In particular the IVC G139.6 + 47.6 
had a relatively large number of background stars 
and few foreground stars. We determine the 
cloud's lower distance limit to be 129 ± 10 pc, 
corresponding to the distance to TYC 4151 1452. 
Our closest detection star is TYC 4152 484, which 
has spectral class F5V, yielding an upper limit of 
2571^ pc. 

Under the assumption that these distances are 
correct and given that every star behind the cloud 
shows a detection, we conclude that the covering 
factor is consistent with unity for this cloud. 

To improve the distance determination, we 
would require more foreground stars that arc more 
distant, as well as background stars that are closer 
than what is currently in our sample. We have 
exhausted the Tycho and USNO catalogs of ap- 
propriate stars. It is possible (though unlikely) 
that more candidates could be found by making 
a photometric survey of the region. By consid- 
ering stars bluer than B — V = 0.35, the survey 
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would have to be complete to apparent magni- 
tudes of my = 10—11 to sample the appropriate 
distances, assuming no extinction. The Tycho 
Catalogue is 99% complete down to my = 11 
(Hocg ct al. 1997) so it is unlikely that we have 
missed candidate stars with spectral type earlier 
than F. 

The distance estimate might also be improved 
by further study of the inconclusive stellar candi- 
dates. In particular, better modeling of the stellar 
features in the template spectra used for the cross- 
correlations could reduce ambiguity and allow one 
to classify a detection or non-detection with more 
confidence. 
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Table 1 
List of Program Stars. 



Tycho ID V MK Class D 

(mag) (parsecs) 



EW(Di) N H i 
EW(D 2 ) (xl0 19 cm- 2 ) 
(mA) 



Notes 



TYC 4152 484 10.45 F5V 

TYC 4385 152 9.66 AOV 

TYC 4148 239 10.91 A1V 

TYC 4151 22 10.53 AOV 

TYC 4152 370 10.27 F0V 

TYC 4145 1284 10.45 F2(IV?) 

TYC 4149 325 9.8 F2(?) 

TYC 4148 1132 8.58 

TYC 4151 1452 6.40 



2571J& 1 
6481"? 
875± 4 7 f 
968±|g| 
342±^ 6 

r 10+421 
O1 ^-207 

229 - 3000 
165 ± 25 
129 ± 10 



150 ± 20 
130 ± 20 
210 ±20 
200 ± 20 
310 ±40 
370 ± 40 
240 ± 40 
270 ± 40 
380 ± 50 
380 ± 50 
330 ± 70 
320 ± 70 
260 ± 50 
250 ±50 



10 detection 

26 detection 

57 detection 

28 detection 

26 possible detection 

37 possible detection 

14 possible detection 

47 possible non-detection 



19 



non-detection 



Note. — EW(Di) and EW(D 2 ) arc the equivalent widths measured using the Na I Di and D 2 lines, 
respectively. 
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Fig. 1. — Positions of our target stars with respect 
to the IVC cloud G139.6 + 47.6. A circle around 
a cross denotes a foreground star whereas a box 
around a cross denotes a background star. Both 
boxes and circles are drawn with an angular size 
inversely proportional to their nominal distances. 
The background gray-scale is from the IRAS lOO^i 
survey. 
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Fig. 2. — H I emission along the line of sight to the 
four stars with detections and the one star with no 
detection (TYC 4151 1452). The velocities of the 
detected absorption are plotted as error bars with 
line styles that match their stars' profiles. 
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Fig. 3. — Intensity of the sky at the DDO, av- 
eraged over one month. Note the self-absorption 
due to the high pressure sodium lights as well as 
the lines due to the low pressure lights. 
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Fig. 4. — Sample spectrum of the telluric standard 
HD177724. 
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Fig. 5. — Na I spectra and cross-correlations 
for the four stars with detections. Two cross- 
correlations are plotted for each star: the sodium 
region, where only the Na I lines are used (solid 
line) and the stellar region, where all lines but the 
Na I are used (dotted line). The cloud velocities 
are over-plotted on the cross-correlations as error 
bars. 
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Fig. 6. — Expected Na I absorption in the spec- 
trum of TYC 4151 1452 if it were behind the 
cloud. The top-most spectrum is the observed 
spectrum. Each subsequent spectrum is labeled 
with the EW(Na I)/Nhi used to predict the equiv- 
alent width one would observe, which is shown in 
parentheses. 
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